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Abstract 

We explore the possibility of distinguishing neutrino mass hierarchies through the neutrino 
signal from dark matter annihilation at neutrino telescopes. We consider a simple extension of the 
standard model where the neutrino masses and mixing angles are obtained via the type-II seesaw 
mechanism as an explicit example. We show that future extensions of IceCube neutrino telescope 
may detect the neutrino signal from DM annihilation at the Galactic Center and inside the Sun, 
and differentiate between the normal and inverted mass hierarchies, in this model. 


1 Introduction 


Although various lines of evidence support the existence of dark matter (DM) in the universe, its 
identity remains as a mystery. Among various proposed candidates, the weakly interacting massive 
particles (WIMPs) are promising ones that typically arises in many extensions of the standard model 
(SM) [1,2]. Major direct, indirect, and collider experimental searches are currently underway to detect 
the particle nature of DM and determine its properties. 

Just like DM, the origin of neutrino mass and mixing defines a very interesting area of investigation 
beyond the SM (BSM) [3,4]. These are the two encouraging avenues of new physics and a large number 
of BSM scenarios have already been proposed to address them. It will be even more interesting to 
investigate models where the two sectors are connected. In fact, such a connection arises naturally 
in a class of models where DM is tied to the neutrino sector. For example, in type-II seesaw plus 
a singlet scalar scenario [5], it has been shown that the neutrino mass hierarchy can influence the 
DM annihilation rate to charged leptons as well as neutrinos and subsequently provide a possible 
explanation of the positron fraction excess observed at the AMS-02 experiment [6]. Similar studies 
have been done in the B — L extension of the minimal supersymmetric standard model (MSSM) where 
it is shown that the annihilation of the right-handed sneutrino DM produces neutrino final states with 
a distinct feature [7]. There have been studies on the possibility of discovering the B — L sneutrino 
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DM at the LHC [8-10]. The connection between the DM and neutrinos in such scenarios can also be 
probed through DM direct detection [11] and some indirect detection experiments [12], Therefore, a 
combination of direct, indirect and collider signatures could be useful to explore this connection (for a 
review, see [13]). 

Our focus here is on the DM indirect detection signals as a means of testing the connection between 
the DM and neutrinos. Neutrino telescopes like IceCube are able to examine these models through 
the signal arising from the DM annihilation into neutrinos at the Galactic Center [14] and inside the 
Sun [15]. Unlike the photon flux, the neutrino flux arising from the DM annihilation at the Galactic 
Center has less astrophysical background uncertainty, which allows us to probe the exact nature of 
the DM more accurately. However, the model prediction for the photon fluxes coming from the DM 
annihilation at the Galactic Center should be consistent with the Fermi-LAT [16] data on the gamma- 
ray flux. Therefore, it will be interesting to combine the results from IceCube and Fermi-LAT to study 
models where the DM and the neutrino sectors are interconnected. 

In this paper, we study these issues within the above mentioned extension of the SM [5,6]. In this 
model, the light neutrino masses and mixing angles are generated using the well-known type-H see¬ 
saw mechanism [17-21] that introduces SU{‘2)l triplet scalars. The neutral component of the triplet 
acquires a non-zero vacuum expectation value (vev) that generates tiny neutrino masses by breaking 
lepton number by two units and mixing among different neutrino flavors. Among the SM particles, 
these triplet scalars only couple to gauge bosons through gauge couplings and to leptons via Yukawa 
couplings. These Yukawa couplings are related to the triplet vev, which in turn determines different 
decay modes of these triplet scalars. It has been observed [22,23] that the triplets will dominantly 
decay to a pair of leptons if the vev is less than 0.1 MeV, otherwise the gauge boson final state will 
become dominant. In order to accommodate a stable DM candidate, the model is augmented by a SM 
singlet scalar with a Z 2 symmetry [5]. In the minimal case, the DM only couples to the Higgs and 
to the triplet scalars. For a vev smaller than 0.1 MeV, the annihilation of the DM produces a pair 
of triplets that will further decay to SM leptons (including neutrinos) with almost 100% branching 
fraction. Therefore, the flavor composition of the final states will depend upon the neutrino mass 
hierarchy that sets the Yukawa couplings. We exploit this feature and show how different neutrino 
mass hierarchies can be distinguished by using the muon signal arising from conversion of muon neu¬ 
trinos from DM annihilation at the Galactic Center and inside the Sun at IceCube. We note that the 
photon signal from the Galactic Center does not discriminate amongst different hierarchies because 
the difference will not be significant. 

The paper is organized as follows. In Sec. 2, we discuss the model in brief. In Sec. 3, we present the 
analysis of the neutrino signal from DM annihilation at the Galactic Center and inside the Sun for 
normal and inverted hierarchies. We also calculate the corresponding muon spectra at the detector 
and compare them with the background. Finally, in Sec. 4 we present our discussion and conclusion. 


2 Type-II seesaw with a scalar singlet DM candidate 


The minimal version of the type-H seesaw model is an extension of the SM that includes a SU{2)l 
scalar triplet A with hypercharge Y = 2. 


A = 


A+ 

C2 

AO 



( 1 ) 
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In order to accommodate a cold dark matter (CDM) candidate, we introduce a SM singlet real scalar 
D. The stability of the DM candidate is ensured by imposing a Z 2 symmetry under which D is odd 
while all other SM particles and the triplet A are even. 


The scalar potential of the model, including the relevant terms for the DM candidate D, is given by [6] 
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Vdm(^,A,D) = -m^D2 + ADD^ + A<jD2($t$) + AAD2TV(AtA). 


( 2 ) 

( 3 ) 


where is the SM Higgs doublet. The couplings Aj(i = 1 — 5) can be taken as real without any loss of 
generality. Due to spontaneous symmetry breaking, the non-zero vev of the Higgs doublet generates a 
tadpole term through interaction involving the Ag coupling in Eq. (2). This in turn induces a non-zero 
vev (ua) for the neutral component of the triplet thereby breaking lepton number by two units. This 
triplet vev contributes to the SM gauge boson masses at the tree level which leads to a deviation in the 
electroweak p parameter. Hence, ua is constrained by the electroweak precision data that requires the 
p parameter to be close to the SM value of unity [24] which eventually demands f a < 2 GeV [25]. In 
addition to this, current experimental bounds on lepton flavor violating processes put a lower bound 
on UA [26-28] as 


uaMa > 150 eV GeV. 


( 4 ) 


It is to be mentioned here that we assume negligible mixing between the doublet and triplet scalars 
in our analysis ® and with this assumption, the DM potential of Eq. (3) can be expressed in terms of 
physical scalars {h, , A^) as 


VDM(4>,A,iA) — 

XaD^ H++H~- + H+H- + I {Hi + Al) + \vaHo 


( 5 ) 


Here, + X^v‘^ -|- Aa^'a denotes the physical mass of the DM candidate. It is evident 

from Eq. (5) that the DM candidate can annihilate to a pair of Higgs and to a pair of exotic triplet 
scalar particles through the coupling parameters A^ and Aa- In the limit moM > niA; where, rriA 
is the degenerate mass of all triplet scalars, the annihilation cross section of the DM candidate (non- 
relativistic) is expressed as 


{av) 
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( 6 ) 


This should yield the correct relic abundance of the DM particle that lies within the 3 ct limit of current 
Planck bound = 0.1199 ± 0.0027 [30]. Now, the final state of DM annihilation will depend 

upon the branching fraction of the Higgs and the triplet scalar decay. The triplet scalars can couple to 
the SM gauge bosons through gauge coupling and to SM lepton doublets through Yukawa couplings. 


Cy = - ^{YA)ijLjCia2ALj + li.c. 

^Detailed expressions for physical scalar mass eigenstates can be found in [29]. 
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Here, C is the charge conjugation operator and Li = is the SU{2)l lepton doublet with i being 

the three generation indices. Further, these Yukawa couplings can be obtained from the Majorana mass 
matrix of neutrinos that arises due to the non-zero triplet vev v/\. 



( 8 ) 

Ya = — =, 

VA VA 

(9) 


where = diag(mi,m 2 ,m 3 ) and U is the PMNS mixing matrix. 

With the recent global analysis data (see Eq. (18) in the Appendix), and after using Eq. (9), we obtain 
the following structure of the Yukawa coupling matrix for the normal hierarchy (NH) and inverted 
hierarchy (IH) scenarios 


= 


10-2 eV 
VA 


X 


1.08 -0.29i 
-1.55 -h0.09i 
1.23-0.31i 
3.84 -h 0.34i 
1.21 -0.13i 
-l.39- 0.94z 


-1.55-F0.09z 
2.07-F0.26i 
-l.59-0.21z 
1.21 -0.13i 
2.97-0.35i 
1.98 -0.65i 


1.23-0.31Z 
-l.59-0.21z 
2.59-FO.Oz 
-l.39- 0.94z 
1.98 -0.65i 
2.66-FO.Oli 


(normal hierarchy) 

( 10 ) 

(inverted hierarchy). 


To find the Yukawa matrix for the degenerate case, we impose the 95% C.L. limit on the sum of all 
light neutrino masses < 0-23 eV from Planck [30], which uses the Baryon Acoustic Oscillation 

(BAO) and Cosmic Microwave Background (CMB) data. By choosing mi ~ m 2 ~ m 3 = 0.07 eV, we 
find 




10-2 eV 

- X 

VA 


6.80 - 0.06i 
-0.13 - 0.04z 
0. - 1.65z 


-0.13-0.04z 
6.91 - 0.03z 
0. - l.lOz 


0. - 1.65z 
0. - l.lOz 
6.71-FO.lOz 


(degenerate case). ( 11 ) 


However, as we will see later, the degenerate scenario does not offer any significant information in our 
study. 

The triplet scalar decay to leptonic final states, —)• —>• , Hq/Aq —)• uv with almost 

100% branching ratio when va < 0.1 MeV. On the other hand, for larger z;a(> 0.1 MeV), gauge 
boson final states become dominant with the decay modes —>■ W~^W^ , H'^ —)• W~^Z , Hq/Aq —)• 

ZZ/W^W~ For almost same order of and Aa within the range for producing correct relic 

density, the DM can annihilate to three different final states (depending upon the branching ratios i.e 
on VA for triplets) : 


(i) neutrino and charged lepton final states. 

(ii) , ZZ finals states mainly from triplet scalar decays; WW*, ZZ* finals states arise from 
the decay of the SM Higgs boson. 

(iii) h , r final states from the decay of the SM Higgs boson. 


We choose two extreme values of va, namely, 1 eV and 1 GeV, respectively to incorporate the effect 
of triplet scalars decay to only leptonic final state or only to non-leptonic final states. Our study is 
mainly focused on the study of neutrino flux coming from DM annihilation which can be obtained from 
all the above three final states. However, for triplet scalars decaying into leptonic final states includes 

^Throughout this analysis we assume degenerate triplet scalars. 
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the contribution of Yukawa couplings Ya which is a function of neutrino mass hierarchy and we expect 
to see the difference in the neutrino flux between the NH and IH scenarios. In the following sections, 
we will explicitly show how our analysis depends upon the neutrino mass hierarchy. 

In this regard, it is to be noted that the triplet scalars can produce interesting signals at the colliders. 
The strongest limit on the scalar masses comes from the current searches at the LHC for the signature 
of doubly charged scalar where they can be pair produced via Drell-Yan and Vector Boson fusion and 
then decay to a pair of same-sign leptons. The experimental lower bound on m}j±± has been set by the 
ATLAS experiment [31] from a pair of isolated lepton search for —)• e^e^) and —?■ 

decay modes at the center of mass energy 8 TeV. The 95% C.L. lower limit on the doubly charged 
Higgs mass for the same-sign isolated muons hnal state exclude mass range below 516 GeV. However, 
the experimental lower limit is based on the assumption that the doubly charged scalar decay to the 
dimoun channel with 100% branching ratio which is not the case in our scenario. In this model, the 
doubly charged scalar decay to dimuon channel with atmost 30% branching ratio and thus the lower 
limit on mu±± can be relaxed. Following Fig. 5 of Reference [31], we set the degenerate mass of the 
triplets {m^) as 400 GeV throughout our analysis. 


3 Analysis 


In this section, we demonstrate our hndings on the DM annihilation at the Galactic Genter and inside 
the Sun and the possibilities of using the respective neutrino signals to distinguish the NH and IH 
scenarios. For the purpose of our analysis, we choose to work with four benchmark points, shown 
in Table 1, by hxing for two DM masses. We choose the DM mass to be greater than 400 GeV 
to get on-shell triplet scalars in the pair production processes of DM annihilation. Our choice of ua 
has already been justihed in the previous section. We choose Aa and couplings in the model (see 
Eq. (3)) such that the correct DM abundance is obtained via thermal freeze-out with nominal value 
for thermally averaged annihilation cross-section (ciann^) — 3 x 10“^® cm^/s. A$ also enters into the 
direct detection cross section {(tsi) through the f-channel exchange of the SM Higgs. Values of A^ 
that satisfy thermal relic abundance yield a direct detection cross section well below the current limits 
from LUX experiment [32]. The choice of A a and A$ and the corresponding relic density and direct 
detection cross section are also listed in Table 1 . The correct DM relic abundance can also be produced 
via non-thermal mechanisms [33-35] in which case a larger DM annihilation cross section is allowed. 
We will discuss the implications of a larger annihilaiton cross section at the end of this section. 

As mentioned earlier, the DM annihilation will give rise to neutrino fluxes of different flavors that 
will help in distinguishing the structure of Yukawa coupling and, hence, the neutrino mass hierarchy. 
The IceGube detects neutrinos by recording the Gherenkov light from relativistic charged particles 
in its volume. Muon neutrinos {v^) produce muon (^) tracks via the charged current interactions in 
the detector. On the other hand, electron neutrinos {ve) and tan neutrinos {vr) result in hadronic 
and electromagnetic cascade events in the ice. Since the cascades are localized, they carry no direc¬ 
tional information, and hence are not good enough for performing a meaningful DM search over the 
background. For this reason, we focus on fluxes that arrive on the Earth. We also estimate the 
photon flux for all scenarios in our model. Gosmic ray showers create a muon background that can be 
controlled by selecting only the upward going events since muons are stopped in the Earth. This limits 
the observation of DM signal to the time when the source is below the horizon. With atmospheric 
muons thus eliminated, the most significant contribution to the remaining background comes from 
atmospheric neutrinos [36]. In addition to this, a portion of the detector may be used as a veto to 
observe the contained muon events with a conversion vertex inside the instrumented volume, as in the 
case of DeepGore array in IceGube. The veto procedure virtually eliminates the contribution to the 
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Benchmarks 

VA 

m.DM(GeV) 

mA(GeV) 

Aa 

'^0 

relic density 

osi (pb) 

BPl 

1 eV 

500 

400 

0.055 

0.04 

0.117 

2.25 X lO-i^' 

BP2 

1 eV 

700 

400 

0.075 

0.05 

0.122 

1.79 X lO-i^' 

BP3 

1 GeV 

500 

400 

0.055 

0.04 

0.117 

2.25 X lO-i^' 

BP4 

1 GeV 

700 

400 

0.075 

0.05 

0.122 

1.79 X 10“^^ 


Table 1: Benchmark Points (BPs) and the corresponding DM relic density and direct detection cross- 
sections. 


background from atmospheric muons by selecting only contained vertices. This increases the potential 
observation time to the full year when the source is both above and below the horizon. 

We use both contained and through-going muons in our analysis. For contained muon tracks, the vertex 
at which is converted to a muon, is within the detector volume. Through-going muons represent 
those events that go through a surface inside the detector but may have been produced outside it. To 
acquire the muon spectra, we use three simple methods: 

1. WimpSim [37] automatically provides muons converted from solar neutrinos. 

2. DarkSUSY [38] gives muons from atmospheric neutrinos (which we average to make an isotropic 
approximation). 

3. Finally, GENIE [39] provides muons converted from neutrinos coming from the Galactic Genter. 

GENIE provides the spectrum of contained muons produced from arriving at the detector. To find 
the total number of muon events, we use the total deep-inelastic charged-current cross sections from 
reference [40] to find the average and P^-to-/i cross sections. Modeling the Galactic Genter as 

an isotropic circle with radius of 5° and ignoring detector angular resolution, we obtain a 5° optimal 
cut, whereas the solar muons have a 2° optimal cut. To find the through-going muons, we propagate 
the muons using the parameters of Table (4) of reference [41]. 

In order to generate the energy spectra of the SM particles that are produced at the annihilation point, 
we first create the model file for GalcHEP using FeynRules 2.0 [42] and then use GalcHEP 3.6.23 [43] 
and Pythia-6 [44] to get the spectra. These spectra are valid both at the Galactic Genter and inside 
the Sun. In the following we present our result for these two distinct DM annihilation points. 


3.1 Signals from DM annihilation in the Galactic Center 

To calculate the fluxes near the Galactic Genter, we run the indirect detection module of mi- 
crOMEGAsv4.1.8 [45] and measure the fluxes for angle of sight 0° < T < 5°. We use the Navarro- 
Frenk-White (NFW) [46] profile p{r) = pQ{r/rs)~^/{'^ + r/rgf' with = 20 kpc, po = 0.4GeV/cm^. It 
is noteworthy that if the annihilation of our proposed DM candidate happens near the Galactic Genter, 
then the final state neutrinos will encounter oscillations on their way to the Earth. In the Appendix, 
we show the flavor dependent probability of neutrino flux that originates at the Galactic Genter and 
reaches the detector surface on the Earth. Hence, the final neutrino fluxes will be given by Eq.(16-19b). 
We have checked that varying all neutrino oscillation parameters within the 3(T allowed range about 
their best fit central value results in 0(3 — 4%) change in the neutrino flux from DM annihilation at 
the Galactic Genter. This is much smaller than the difference due to different mass hierarchies. 

Following are the interesting features that are seen in Figs. 1, 2 and 3. 
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• Fig. 1 displays the fluxes for NH and IH for BPl and BP2 which is for ua = 1 eV. As 
previously argued, at such low vev, the triplet scalars generate leptonic final state with almost 
100% branching ratio and the DM annihilation to two on-shell triplet scalars mainly contribute to 
the annihilation cross section, hence to the flux. Now, it is to be observed that for moM = 500 
GeV, the flux rises at around 100 GeV and then falls near 400 GeV and similarly, the rise and fall 
occur at 60 GeV and 640 GeV for tudm = 700 GeV. This is not surprising and can be understood 
by the kinematics of the DM annihilation. The triplets are produced on-shell with some boost 
that comes from the mass difference between the DM and the triplets. Thus the two-body final 
states render a box like feature with end points mDM(l 4= /3)/2 where j3 = y^l — rn?^/ 

• We also see from Fig. 1 that there exists a significant difference between the NH and the IH 
scenarios. This is due to the fact that in the NH case, we get more taus (r) from the triplet 
decays which produces and the Vr further gets converted into For both the DM mass, 
the neutrino flux can be easily a factor of two higher in the NH scenario than the IH case. The 
neutrino flux due to degenerate case will be in between the flux generated for NH and IH cases. 
This is simply because in the degenerate case, we get less rs from the triplet decays compared to 
the NH case but more than the IH case. This is true for all our discussions to follow. Therefore, 
the degenerate case will not provide any new insight and so we restrain ourselves from doing any 
further remark on this case. 



£,.,(GeV) 



Figure 1: Spectra of from DM annihilation at the Galactic Genter for BPl (left) with ttidm = 
500 GeV and BP2 (right) with ruDM = 700 GeV. 


• In Fig. 2, we present the diffused photon fluxes for NH and IH cases. In our model, the photon 
flux arises from external charged lepton legs, final state radiation, and secondary decay of charged 
leptons that are directly produced from triplet decay. For both the benchmark points, we observe 
that the photon flux does not distinguish between NH and IH cases 

• With the increase in the branching ratio of the triplet scalars to leptons reduces and the decay 
modes with gauge boson final states open up and at sufficiently large ua, BR(77^^ —)• W^W^} ~ 
100%, dominating over the leptonic final states. Hence, for ua = 1 GeV, the dependence on 
neutrino mass hierarchy gets washed away. Moreover, the contribution to neutrino fluxes coming 

®The DM annihilation can also produce photon lines via charged scalar loop processes DM -|- DM —>■ 77 and DM -|- 

m2 

DM —>■ Z'y with respective energies tudm and tudmII -However, the line signal is highly suppressed relative to 

the diffused photon signal in our model, which is in agreement with the fact that that Fermi-LAT has not observed any 
photon line at such energies. 
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Figure 2: Photon spectra from DM annihilation at the Galactic Center for BPl (left) with ttidm = 
500 GeV and BP2 (right) with mDM = 700 GeV. 


from DM annihilating to Higgs will become more pronounced in this case. In Fig. 3, we show 
neutrino fluxes for BP3 and BP4 and evidently these do not distinguish between NH and IH 
cases. However, the shapes of neutrino spectra are different compared to the previous ua = 1 eV 
cases as shown in the Fig. 1. As expected, the flux is higher in the previous case because the 
leptonic final states originate from direct triplet decay. The photon flux, however, does not draw 
any discrimination be it between nA=l eV and 1 GeV cases or between NH and IH scenarios. 



£v,(GeV) 



EyiGeV) 


Figure 3: Spectra of (left) and Photon (right) from DM annihilation at the Galactic Center for BP3 
and BP4. 


3.2 Signals from DM annihilations inside the Sun 

DM particles that pass through the Sun lose energy due to the DM-nucleon scattering and get gravi¬ 
tationally trapped [47-53]. The total rate of DM annihilations Fann inside the Sun follows [1] 

Hann = ^tanh^{VCAt). ( 12 ) 

Here C is the capture rate of DM particles by the Sun, which is related to crgi in our model, and A is 
related to the DM annihilation cross section dann (for details, see [1]). DM capture and annihilation 
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reach equilibrium when t > Teq = {VCA)~^^ in which case Tann ~ C/2. For values of cisi given in 
Table 1, and (cann'i^) = 3 x 10“^® cm^ s“^, the equilibrium condition is achieved for the age of the Sun. 
We note that as long as equilibrium is established, the total rate of DM annihilations is controlled 
by the direct detection cross section cjsi. Therefore, increasing (Tann (which is allowed in scenarios of 
non-thermal DM) and keeping usi hxed does not enhance the signal from DM annihilation inside the 
Sun. It only leads to a faster establishment of equilibrium between DM capture and annihilation. 

DM annihilation hnal states that are most relevant for producing neutrinos inside the Sun are prompt 
neutrinos from /A^ and /H~ decays, W, Z and r particles from /H~ and /H decays. 

Other hnal states (like e, and lighter quarks) lose their energy and stop almost immediately, which re¬ 
sults in production of neutrinos at energies below detection threshold. The charged current interactions 
inside the Sun convert Ug, Ur to e, /r, r respectively. Using the charged current neutrino-nucleon 
cross section [40], we hnd that neutrino absorption becomes important (i.e.. Tabs < Rci with Rq being 
the core radius of the Sun) at energies > 300 GeV. For I'e, the havor and mass eigenstates are 
the same deep inside the Sun, and hence absorption suppresses the hux of at energies above 300 
GeV. On the other hand, ly-j- absorption produces r particles that decay quickly before losing too much 
energy because of very short lifetime of r. This decay produces Vr at a lower energy, and hence this 
“regeneration” effect populates the spectrum at energies well below the DM mass. Since and 
undergo oscillations inside the Sun that is dominantly set by the atmospheric mass splitting 
we have Lose ot long as L^bs ^ Losc/4, oscillations mix and efficiently and z^^ hnal 

states also feel the regeneration effect. When Labs drops below Losc/4, which happens at ~ 500 
GeV, oscillations cease to be effective and gets absorbed similar to Ve- As a consequence, only the 
Ur hnal state retains a signiheant regeneration signature at energies above 500 GeV.® 

For neutrinos of sufficiently low energies, vacuum havor oscillations between the Sun and the Earth 
are averaged over half a year when the Sun is below the horizon at the south pole. In particular, below 
100 GeV, the oscillation length set by the solar mass splitting is less than approximately 3 

million kilometer change in the Earth-Sun distance over half a year. The situation changes at energies 
above 100 GeV where solar neutrino oscillations are not averaged out anymore. For monochromatic 
neutrinos of a single havor, this signihcantly affects oscillations between Ug and u^jur-, and hence the 
u^ spectrum at the detector, at high energies. However, this is not an important effect in our model 
since DM annihilation produces different neutrino havors with a continuous energy spectrum in this 
case. 

We use DarkSUSY 5.1.1 [38] to simulate production of neutrinos in the Sun, their propagation to 
the South Pole, and the interaction of Vy, with ice at many different energies for all the SM channels. 
DarkSUSY does this by interpolating an older WimpSim [37] simulation. Here we run WimpSim 3.05 
directly for the prompt neutrino channels , sample energies spaced at 10 GeV, and use oscillation 
parameters found in [54]. 

We then combine the energy spectra of all the channels at the production point with how they prop¬ 
agate, giving us the hnal neutrino spectra at the detector.^® It should be noted that variation of all 
neutrino oscillation parameters within the 3(T allowed range results in 0(3 — 4%) change in the hnal 
neutrino spectra from the Sun, similar to the Galactic Genter case, which is insigniheant. 

In Fig. 4 we show the spectra of Uy at the detector for ua = 1 eV (BPl, BP2) and for ua = 1 GeV 
(BP3) cases. It is seen that at energies below 400 GeV the NH scenario results in a larger number of Uy 

^Neutrinos also have neutral current interactions with matter inside the Sun that results in energy loss of the neutrinos 
from all flavors further shifting their spectra toward lower energies. However, the cross section for neutral current 
interactions is a factor of 3 smaller than that for charged current interactions, which makes them subdominant. 

^°Since i/-nucleon cross section is about two times smaller than the zz-nucleon cross section [55,56], the number of 
absorbed v inside the Sun is smaller. This results in a larger v flux at the detector from DM annihilation inside the Sun. 
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and eventually muons than the IH scenario. The reason is similar to the Galactic Center case and can 
be explained by the triplet scalar decay patterns. For ua = 1 eV, as mentioned above, triplet decays 
mainly produce Vr, and r, ^ final states in the NH scenario. In the IH scenario, on the other hand, 
triplet decays mainly produce and e final states. As far as charged lepton final states are concerned, 
only r decays before losing a significant fraction of energy and produces neutrinos. Regarding triplet 
decays to neutrinos, Vr (and to a lesser extent v^) is the most relevant flavor that survive at energies 
above 300 GeV due to regeneration effect. In consequence, for = 1 eV, a larger number of 
arrive at the detector in the NH scenario. However, the neutrino spectra in the two scenarios do not 
differ when = 1 GeV as in earlier case. Gomparing the neutrino signal from DM annihilation in 




A^GeV) 


Figure 4: Spectra of from DM annihilation inside the Sun at the detector for tudm = 
500 GeV (700 GeV) in the left (right) panel. 

the Sun (Fig. 4) and at the Galactic Genter (Fig. 1), we notice some differences. First, the kinematic 
cuts in Fig. 1 are not prominent in Fig. 4. This is due to the fact that absorption/regeneration and 
down scattering inside the Sun erases sharp features in the spectra. More precisely, only a little bump 
develops at low energies for moM = 500 GeV due to the boosted triplet scalars, shown in right panel 
of Fig. 4, unlike the box-like structure that can clearly be observed in the Galactic Genter neutrino 
flux. At higher energies, instead of a bump only the wiggles can be noticed which are the effect of 
neutrino oscillation. Furthermore, the difference between the spectra in NH and IH scenarios in Fig. 4 
is smaller than that in Fig. 1. This is because stopping of muons and partial absoprtion of inside 
the Sun decreases the overall flux of arriving at the detector. Nevertheless, the NH scenario yields 
a larger flux than the IH scenario when ua = 1 eV. 

We emphasize that the neutrino signal from DM annihilation inside the Sun is complementary to that 
from the Galactic DM annihilation as they are set by different cross sections (ctsi and (Tann respectively). 
We see from Figs. 1 and 4 that for the model parameters given in Table 1 the signal from the Sun is 
about an order of magnitude stronger than that from the Galactic Genter. As we will discuss later, 
both of these signals can be enhanced further. 


3.3 Muon spectra at the detector 

We now present the muon spectra obtained from the conversion of at the detector, which is the 
observed signal at neutrino telescopes. For simulation purpose, we assume a detector that has the 
same capability as the IceGube DeepGore array. In Fig. 5, we show the spectra of contained and 
through-going muons in the detector from DM annihilation at the Galactic Genter. In this case, the 
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Figure 5: Spectra of muons from DM annihilation at the Galactic Center for m-oM = 
500 GeV (700 GeV) in the left (right) panel. Upper and lower panels show contained and through-going 
muon spectra with an angular cut of 5° respectively. 


signal comes from a region around the Galactic Center that has 5° angular extension. Hence, in order 
to optimize the signal to background ratio, we have imposed a 5° angular cut on the muons relative 
to the center of the galaxy. 

In Fig. 6, we show the spectra of contained and through-going muons from DM annihilation inside 
the Sun. Since the Sun is a point-like source, the optimal signal to background ratio is obtained for 
an angular cut of 2°, in this case. The peaks in both figures are due to the imposed angular cuts 
that eliminate muons produced from conversion below a certain energy. For through-going muons 
the peak occurs at a lower energy, which can be understood by noting that the measured energy of a 
through-going muon is in general less than the actual energy at the production point. 

As expected, the muon spectra have less features than the neutrino spectra shown in Figs. 1-4. Nev¬ 
ertheless, the difference between ua = 1 eV and ua = 1 GeV cases, as well as the NH and IH scenarios 
in the former case, are clearly visible in the muon spectra arising from both the Galactic Genter and 
the Sun. Similar to the neutrino spectra, substantial deviation among the cases is seen in the signal 
arising from DM annihilation at the Galactic Genter than in the Sun. However, we also note that the 
absolute value of the flux from the Galactic Genter is smaller by an order of magnitude than that from 
the Sun. 
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Figure 6: Spectra of muons from DM annihilation in the Sun for ruDM = 500 GeV (700 GeV) in the 
left (right) panel. Upper and lower panels show contained and through-going muon spectra with an 
angular cut of 2 ° respectively. 


The muon spectra shown in Figs. 5 and 6 are obtained assuming the nominal thermal freeze-out value 
of {(Jaunv)- However, the Fermi-LAT data from Milky Way dwarf spheroidal galaxies [57] allows larger 
values on the DM annihilation cross section for the values of DM mass that we have considered here. 
One can see from Fig. 8 of [57] that for the bb and final states (uann^) can be larger than 

the thermal freeze-out value by a factor of 4 ( 6 ) when w-dm = 500 (700) GeV. In BPl and BP2, 
where the distinction between the NH and IH scenarios is significant, these final states arise from DM 
annihilation to the SM Higgs that is controlled by the coupling A$ of Eq. (3). As explained before, A$ 
also controls dsi in onr model. Hence raising {aannv) by the above factors, due to a larger value of A#, 
also results in an increase in usi compared with that given in Table 1, which is still consistent with the 
LUX bounds [32]. For the final state, {(Jannv) can be larger than the thermal freeze-out value 

by a factor of 10 (17) when tudm = 500 (700) GeV This final state arises due to DM annihilation 
to the triplet Higgs, which is controlled by the coupling Aa of Eq. (3). Increasing Aa results in an 
increase in a^nn, but does not affect usi. 

Therefore, by invoking non-thermal mechanisms for DM production in the early universe, we can obtain 
considerably larger neutrino fluxes from DM annihilation at the Galactic Center and inside the Sun. 

It should be noted that the Fermi-LAT bound cannot be directly applied to our model since we have four-body 
annihilation final states in this case. Moreover, we have a combination of different hnal states instead of 100% of just 
one hnal state. Nevertheless, the Fermi-LAT limit provides a reasonable upper bound on cTann in this case too. 
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The resulting enhancement in the muon events, compared to Figs. 5 and 6, leads to a better prospect for 
detection of the neutrino signal against the atmospheric background. In Fig. 7, we show the spectra of 
contained muons due to DM annihilation at the Galactic Center (left panel) and inside the Sun (right 
panel) for annihilation cross sections that are just below the Fermi-LAT limits, where tudm = 500 
GeV and the background from atmospheric neutrinos is also shown for comparison. Contained muons 
from DM annihilation inside the Sun provide the best detection opportunity. Within energy interval 
100-400 GeV, the total number of contained muons for the NH and IH scenarios is 11 and 9 respectively 
(compared with 95 for the background). The maximum signal to background ratio, 17% and 15% for 
the NH and IH scenarios respectively, occurs at an energy of 255 GeV. For a neutrino telescope with 
the same capability as the IceCube DeepCore array, 3a discovery of NH and IH scenarios takes 8 and 
12 years respectively. To distinguish the different neutrino mass hierarchies, one should go beyond the 
simple number count and perform a careful shape analysis, which is beyond the scope of this paper. 

There are proposals for directly measuring the mass hierarchy by using atmospheric neutrinos in future 
extensions of neutrino telescopes such as PINGU (Precision IceGube Next Generation Upgrade) [58]. It 
is possible to exclude the wrong mass ordering by this method (as well as in neutrino beam experiments) 
at the 3cj level within the next 10-15 years [59]. Our approach, which exploits an interesting connection 
between DM and neutrinos, is complementary to this direct method and makes another case for the 
future neutrino telescopes with much larger volume like IceGube-Gen2 [60]. 




Figure 7: Spectra of contained muons from DM annihilation at the Galactic Genter (left) and inside 
the Sun (right) for moM = 500 GeV. The value of fiann used to obtain the muon flux is just below the 
current bounds from Fermi-LAT. The top line in both hgures is due to the background arising from 
the atmospheric neutrinos. 


4 Summary and conclusions 


In this paper, we have considered an extension of the SM that explains neutrino masses and mixing 
angles via the type-H seesaw mechanism and includes a singlet scalar as a viable DM candidate. The 
DM interacts with the triplet scalar that generates light neutrino masses thereby linking the two 
sectors. We have studied the neutrino signal from DM annihilation at the Galactic Genter and inside 
the Sun in this model that could be detected at IceGube. Our main results are summarized as follows: 

• In the region of the parameter space where the triplet scalar dominantly decays to leptonic final 
states, the flux of from DM annihilation depends upon the neutrino mass hierarchy (normal. 
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inverted or degenerate). The photon flux, on the other hand, is practically the same for different 
hierarchies. 

• The difference in the flux of at the detector for NH and IH is visible. The NH scenario produces 
a larger flux because of its mass pattern. This holds for both signals from DM annihilation at 
the Galactic Center and inside the Sun. 

• The difference is more significant in the neutrino signal arising from DM annihilation at the 
Galactic Center than that arises from DM annihilation inside the Sun. This is mainly because 
the interaction of neutrinos with matter inside the Sun results in a moderate attenuation of the 
final flux at the detector. 

• The same reason holds for the contained and through-going muons that are produced from the 
conversion of at the detector. The muon spectra have less features than the neutrino spectra. 
Nevertheless, these muons can be detected over the atmospheric background. 

• For annihilation cross sections just below the Fermi-LAT limits, future extensions of IceCube are 
able to discover the neutrino signal and distinguish the different neutrino mass hierarchies with 
multiyear data. 

The LHC in tandem with direct detection experiments, the IceCube and Fermi-LAT will be able to 
probe this minimal extension of the SM that utilizes a new Higgs sector. Finally, we would like to 
conclude by pointing out the fact that this analysis can be applied to other models that link DM to 
the neutrino sector. 
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Appendix 


A Oscillation Probability Calculation 


The flavor eigenstate of neutrinos can be related to the mass eigenstate by the 3x3 Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) mixing matrix. 
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The neutrinos from the Galactic Center encounter flavor oscillation while reaching the Earth. However, 
since it traverses very large distance with respect to its oscillation length, hence the oscillation term is 
averaged out. The probability of one neutrino flavor to oscillate into other becomes simply. 
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The PMNS matrix is usually parametrized in terms of the three mixing angles 0i2;023;0i3 j and one 
Dirac (5) and two Majorana (Q!i;a 2 ) CP phases: 


U = 


C12C13 'S12C13 si3e“‘" 

—'SI 2 C 23 — Cl2S23'Sl3e*'^ C 12 C 23 — Sl2S23'Sl3e*'^ 'S23C13 

•Sl2'S23 — Cl2C23Sl3e*‘^ “C12S23 “ 'Sl2C23Sl3e*^ C23C13 


X diag(e*“C2^e*“2/^l) (15) 


The oscillated flux of the neutrinos at the detector at earth. 


/ 4>, \ / Pn Pi2 Pi3 \ f 

= P21 P22 P23 (16) 

\ P31 P32 P33 ) \ ) 

From Eq. (14b) and Eq. (14c), 

Pll = \Ue,\''+\Ue,\^ + \Ue.f (17a) 

Pl2 = \Ue,\^\U^,\^ + \Ue,\^\U^,\^ + \Ue,\^\U^,\^ = P2l (17b) 

Pl3 = \UeA^\U^A^ + \Ue^\^\Ur^\^ + \Ue,\‘^\Ur,\^ = P^i (17c) 

P 22 = \U^X + \U^,\^ + \U^,\^ (17d) 

P 23 = \U^,AUr,\^ + \U^,\^\Ur,\^ + \U^,,\^\Ur,\^ = P ^2 (17e) 

7^33 = \UrA^ + \Ur,\^ + \Ur,\^ (17f) 


We use the latest global analysis of the 3-neutrino oscillation data [54], 


= 7.60 X 10“® eV^ , 
Amii = 2.48 X 10"^ eV^ (NH), 
h = 1 . 4 l 7 r (NH), 
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(IH). 


j 

f 0.5380 

0.1979 

0.2640 

Normal Hierarchy : 

P = 

0.1979 

0.4240 

0.3779 


1 

0.2640 
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j 
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Inverted Hierarchy : 
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